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Abstract. We develop a framework for a unified treatment of well-posedness for the Stefan prob- 
lem with and without surface tension. In the absence of surface tension, we provide new estimates 
for the regularity of the moving free-surface. We construct solutions as a limit of a carefully chosen 
sequence of approximate solutions to our so-called K-problem, in which moving surface is regular- 
ized and the boundary condition is modified. We conclude by proving that solutions of the Stefan 
problem with positive surface tension a converge to solutions of the classical Stefan problem as 

1. Introduction 

1.1. The problem formulation. We consider the local weh-posedness and boundary regularity 
of solutions to the classical one-phase Stefan problem, describing the evolving phase boundary of a 
freezing liquid. The temperature of the liquid p{t,x) and the a priori unknown moving boundary 
r(t) must satisfy the following system of equations: 

(1.1a) pt^ 
(1.1b) 
(1.1c) 

(i.id) p(o,.)=po,r(o) = ro. 

Here, n{t) is an evolving open subset of R'' with d>2, with r(i) denoting the moving boundary 
(which may be a connected subset of dfl{t) if a part of the boundary is fixed). 

Equation (|l.lap expresses the fact that heat heat diffuses in the bulk i}{t), while the boundary 
condition (|l.lb|) states that the heat flux across the boundary governs the boundary evolution; that 
is, dnP = ^p-n is the normal derivative of p on T(t) where n stands for the outward pointing unit 
normal, and VY(t) denotes the speed or the normal velocity of the hypersurface T{t). In the case 
that (7 = 0, (ll.lcl) is termed the classical Stefan condition and problem (jl.ll) is called the classical 
Stefan problem. In this case, freezing of the liquid occurs at a constant temperature p = 0. On the 
other hand, if cr > in (|l.lc|) then the boundary condition is called the Gibbs-Thomson correction to 
the classical Stefan condition, and the system (II. ip is then termed the Stefan problem with surface 
tension, whereby a > is a given coefficient of surface tension and Kr{t) stands for the mean curvature 
of the moving boundary T{t). Finally, we equip the problem with suitable initial conditions (|l.ldl) : 
Po : ^^(0) — >■ K and Fq are the prescribed initial temperature and boundary, respectively. 

Problem (II. ip falls into the category of free boundary problems that require further assumptions 
on the initial data in order to be well-posed. In our context, such a stability condition is called the 
Taylor sign condition and it reads: 

dnPQ>0 onr(O). 

We will come back to this later in the introduction and explain how it naturally emerges from our 
analysis. 

1.2. The reference domain fl and the dimension. We will begin the analysis with motion 
in K^, and then describe the minor modifications which arise for motion in M.^. Also, in order to 
simplify our presentation, we will choose our initial domain to be 
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where T"^ stands for the 1-tonis. This allows us to use one global Cartesian coordinate system 
(rather than a collection of local coordinate charts) , which is ideal for describing new identities that 
provide very natural estimates for the second-fundamental form of the evolving free-boundary r(t). 
All of our results apply to general domains, however, via a finite covering of fl with local coordinate 
charts and a partition-of-unity subordinate to this cover; in particular, the Stefan problem localizes 
to each chart, and effectively reduces to the analysis on r2 = T^x]0,l[. 
At initial time t = 0, the bottom boundary 

T = T^x{x'^ = 0} 

represents the initial reference boundary of the moving domain T{t). The top boundary 90top = 
X {x^ = 1} is fixed in time, and the temperature p satisfies the homogeneous Neumann boundary 
condition on dil: 

(1.2) dnP = on SQtop- 

1.3. Notation. For functions defined on the whole domain il, we generally use the double bar 
notation to designate its various Sobolev norms, and for those defined on F, we use the single bar 
notation. In other words, for any s > and given functions / : f2 — >■ M, : F — )■ M we set 

\\f\\s-=^\\f\\H^{n)] \ip\s:=M\H'{r)- 

In particular, when s is not an integer, the corresponding fractional Sobolev space is defined by 
interpolation in a standard way. If /: [0,r] x O— >-]R, (p: [0,T]xF— >-R are given time-dependent 
functions, then for any 1 < p < oo we set 

ll/llLfHI := ll/l|LP([0,T];H»(n))||; I^PIl^H^ ■= II VlUp([0,T];H»(r)) • 

If i = l,2 then f,i:=dxif is the partial derivative of / with respect to coordinate. Similarly, 
f,ij d^idxi f and so on. When differentiating with respect to the time variable t, we often drop the 
comma, i.e. ft = f,t = dtf- In order to denote the tangential derivatives we introduce the following 
notation: 

df:=f,i, aV:=a^i/. 

We use C to denote a universal constant that may vary from line to line. In numerous estimates 
the sign < is used: by definition X <Y if and only if there exists a universal constant C such that 
X < CY. We use P to denote a generic real polynomial with positive coefficients that can similarly 
vary from line to line. We always sum over the repeated indices. 

1.4. The Arbitrary Lagrangian-Eulerian (ALE) formulation. In order to obtain a priori 

estimates, and to facilitate the construction of solutions, we transform the Stefan problem to an 
equivalent problem on a fixed domain. In many situations, the use of Lagrangian coordinates is 
most natural for this purpose, but for the Stefan problem, this is not the case. Rather, we make use 
of the so-called Arbitrary Lagrangian-Eulerian (ALE) coordinates, which interpolate between the 
Lagrangian and Eulerian representations. 
To this end, we represent 

F(t):=*(i,F) 

as a graph over F. We let h{t,-) : F ^ M denote the height function of the graph over F and for a;' e , 
we define the map 

^{t,x') = {x' Mt.x')) 

so that the moving surface T(t) is parameterized as follows: F(f) = {xGR^| ■s.= {x' ,h{t,x'))} where 
h is the unknown. Note that the outward unit normal n{t,x') in these coordinates, by abuse of 
notation reads: 

. ., (dh.-l) 
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Assuming that h is sufficiently regular and remains a graph, we can define a difFeomorphism 5* : 
f2— >-r2(t) as a (nonlinear) elliptic extension of the boundary diffeomorphism {x',h), by solving the 
following nonlinear equation: 

(1.3a) J(t):^detD^'^\n{t)\/\n\ in n, 

(1.3b) '^{t,x',0)^{x',h{t,x')) x'er, 

(1.3c) * = Id on OQtop- 

Hereby \U\ stands for the Lebesgue volume of the set U. Note that J(0) = 1 and for short times, 
J{t) will remain close to 1. The existence and regularity of solutions to (jl.3p is guaranteed by a 
theorem of Dacorogna and Moser [3D], with the important addition that for s>0.5 the following 
elliptic estimate holds: 

(1.4) ||*lk=(a)<q|*lk=-o..(r). 

In fact, Dacorogna and Moser prove the result in Holder spaces, but the argument works in the 
same fashion in Sobolev spaces of integer order, and in Sobolev spaces of fractional order by linear 
interpolation. This inequality (jl.4p provides optimal control for derivatives of ^I^ in 51. We shall 
often use the vector parallel to n given by 

2,_{dh,-l) 



(1.5) Ai:= J , 
where we recaU J{t)^\n{t)\/\D.\ by (fL3a| . 

1.4.1. ALE change of variables. We introduce 

(1.6) g:=po^', A:=[D^Y^, w.^^t, t;:=-Vpo^'. 
Note that the relation v — —Wpo'i' transforms into: 

(1.7) v' + A';q,k = in n. 

Let us express pto'i' in terms of q,v,w. By the chain rule pt — qto'i'^^ +Dqo^i^^ -"^lY^. Since 
*j"^ = -^o\['-i*to4'-i and w = ^t, using ([TT7)) we obtain 

pto^ = qt- A';w'' = qt + v-w. 

Recall that the transformed Laplacian in new coordinates simply reads 

Apo^ = A^q = Al{A'^q,k),j. 

Remark 1 (Differentiation rules). When differentiating the entries of the matrix A=[D"i!]^^, we 
use the following identities: for given i,fce {1,2} we have 

dtA'l^~A';w\,A^; dA^^-A';d^\,Al 

In particular, simple application of the above identities and product rule imply for any given a,b£ N: 
(1.8) 



i+i'>i 



where the term {•, •} is the commutator error. Here the constants ai^ii are some universal constants, 
depending only on m, n, I and I' (where 0<l<m, 0<l' <n). 

1.4.2. Classical Stefan problem in ALE variables. In ALE coordinates, the classical Stefan problem 
i.e. problem ()l.ip with a = 0) on the fixed domain D, is given by the following system of equations: 

1.9a) qt-Aj{A'^q,k),j = -v-w in fl; 

1.9b) v' + A'yq,k = in 

1.9c) q^O on F; 

1.9d) *f •n(t) = -w-n(t) on F; 

1.9e) v-N = on 9fitop; 

1.9f) q(0,.)=g=poovl/; vl/(0,.)=ld(-), 
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where IdiJl— ^fi denotes the identity map ld{x)=x and A^=(0,1) is the unit normal on 9f2top- 
Problem (ll.Qp is the equivalent formulation of the problem (jl.ip . The map 4* evolves according 
to (|1.9dl) which in turn is coupled to the evolution of q via (|1.9ap . Associated to it is the following 
high-order energy functional: 

a+2b<4 a+2&<3 

(1.10) 2 12 1 

b=0 b=0 b=0 6=0 

Remark 2. It is helpful to keep in mind that the above spaces are simply space-time parabolic Sobolev 
spaces, where the formal derivative count works as follows: one time derivative corresponds to two 
spatial derivatives. For instance, q£ H"^, qt£L^H'^, and qtt&L'^L'^. 

All the time integrals in the L^-norms above are over the time- interval [0,t]. In order to obtain a 
locally well-posed problem for arbitrarily large initial data, we must impose a sign condition on the 
normal derivative of p on Tq at time t = 0: 

(1.11) d„PQ>0 onr(O), 
which expressed in terms of q simply reads 

(1.12) -{qo),2\t=o>OonT. 
The condition (|1.12p ensures that on [0,T] 

(1.13) A(t):= inf (-g,2)(t,x',0)>0 

if T > is taken sufficiently small. Note that —q,2 — ctv-n on F with v defined by (II. 6p and a = 
|9^'|/J~1. As mentioned in Section 1 1.1) we shall refer to (|1.12p as the Taylor sign condition in 
analogy to the terminology used in the well-posedness theory in fluid mechanics j 44| 139] . The Taylor 
sign assumption emerges naturally in our analysis as a sign condition that gives positive-definiteness 
of the energy functional. 

We shall give a more detailed description of the various norms in the energy function £ below. 
We define 

(1.14) Sit):={iq,h): £{q,h){t) <^}- 

Our first result is a well-posedness statement for the classical Stefan problem. 

Theorem 1.1. Given initial conditions {qo,ho)€S with qo satisfying the Taylor sign condi- 
tion hl.l^j) . the problem hl.9j) is locally-in-time well-posed, i.e. there is a T>0 such that there 
exists a unique solution {q,h) with the initial data {qo,ho) on the time interval [0,T]. In addition, 
there exists a universal constant C > such that: 

£iq,h)<2£{qo,ho)- 

Remark 3. The definition of our higher-order energy function £ restricted to time t = requires 
an explanation of time-derivates of q and h at t = 0. Specifically, the values qt\t=o^ Qtt\t=o> ^t|t=o 
and htt\t=o o,'"'^ defined via space-derivatives using equations il.9a\) and \1.9d\) . To ensure that 
the solution is continuously differentiable with respect to t at t — we must impose corresponding 
compatibility conditions on the initial data. In particular, restricting lll.9a\) to T and evaluating at 
time t — 0, we arrive at 

(1.15) (go),22+(go),2-0 on r. 

It is an easy exercise to see that the set of initial data satisfying m.l5\) . q — on T and dnq — on 
diltop is non-empty. For example, notice that the function r(x^ ,x'^) =hi(l-\-x'^) satisfies the first two 
conditions in a neighborhood ofT. By employing a partition-of-unity offl, we ensure that the third 
condition is also satisfied. Higher regularity on the initial data necessitates the use of more compat- 
ibility conditions; they arise by successive time- differentiation of our system of equations. Keeping 
in mind \1.15]) and the fact that \1.9d\) at time t = Q reduces to ht =q,2 on T, we differentiate fl.9a\) 
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with respect to time. Evaluating it at time t — and restricting the obtained equation to the boundary 
T, we obtain the second compatibility condition 

(1.16) -«',2t|t=og,22+2*\2t|t=og,2i+2g,2i = on F. 

Since we imposed the homogeneous Neumann condition on the top boundary dVLtop, we impose two 
further compatibility conditions 



(1-17) '?,2|t=o = 0, g,222|t=o on r2 



top • 



Our second main result establishes the vanishing surface tension hmit. If we replace the boundary 
condition (|1.9c|) with 

(1.18) c ^ 



|a2*|2 (l + |a/l|2)3/2' 

then the problem ([0]) . with (fTTSl) in lieu of (fr9c| . is the ALE version of the Stefan problem with 
surface tension. To state the main result, we define the high-order energy adapted to the presence 
of surface tension: 



(1.19) £'^=riq,h)^£iq,h) + aJ2 |5?/^|^^^5-2. +aJ2 \dth 



b=0 6=0 



For a given T > let 
(1.20) 

C''':={{q,h): qeC\[0,T];C°mnC°{[0,T];C\n)), /iG ci([o,r];C"(r))nc°([o,r];c2(r))} 
with the associated norm 

\\q,h\\ci:2 = II qllc^c^ + 11 alleles + l^lc«c^ + Icicles- 

Definition 1 (Well prepared data). A family of initial data (q^ ,hQ)^>o is well-prepared if 
it satisfies 1) the compatibility conditions associated to the Stefan problem with surface ten- 
sion U.9a\) . U.18]) . U.9d\) . U.9e]) . 2) the Taylor sign condition U.12\) and finally 3) £(gQ,ft.Q)— >■ 
£{qo,ho) as a^O. 

Theorem 1.2. Let ((/q ,/iq)o.>o &e a given family of well-prepared initial conditions in the sense of 
Definition\^ By [q'^ ,h'^)a>o we denote the family of solutions to the associated Cauchy problem for 
the Stefan problem with surface tension \1.9a\) . \1.18]) . il.9d\) . U.9e\) . There exists a cr -independent 
time T>0 and a constant C depending only on {qQ,hQ) such that 

E^iq" X){T)<C CT>0. 

for all cr>0. As a consequence, the sequence {q'^^h'^) converges to the unique solution {q,h) of the 
classical Stefan problem m.9\) with a — in the 



Remark 4. Note that Theorem also contains a local-in-time well-posedness statement for the 
Stefan problem with surface tension under the additional assumption that the Taylor sign condition 
holds for the initial datum. That condition is, of course, superfluous if one wants to show only 
local-in-time well-posedness in the presence of surface tension. A straightforward adaptation of 
our analysis would yield local-in-time well-posedness for the one-phase Stefan problem with .surface 
tension, without any .sign condition on — (7,2 (0,a;',0). 



Remark 5. From the proof of Theorem ] 1. II!] it becomes apparent that we can prove convergence in any 
functional space that embeds compactly into our energy space. We chose to express this convergence 
in C^'^-norm as it guarantees that the limiting solution is a classical one. 



Remark 6. As in Theorem \l.l\ we require suitable compatibility conditions on the initial datum. 
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1.5. Prior results and the motivation for the current treatment. There is a large amount 
of literature on the classical one-phase Stefan problem. For a comprehensive overview we refer the 
reader to Meirmanov [864 and Visintin [JJ. First weak solutions were defined by Ladyzhen- 
SKAYA, SOLONNiKOV and Uralceva [53]. These weak solutions were analyzed by Friedman, 
KiNDERLEHRER [U, [^S], [15], Cafarelli, Evans [3, 0, wherein the regularity of weak solutions 
was established. Since the problem satisfies a maximum principle, it is ideally suited to the so- 
called viscosity solutions approach. Existence and regularity of viscosity solutions was established 
by Athanasopoulos, Caffarelli, and Salsa in [4], [5]. For an exhaustive overview and intro- 
duction to the regularity theory of such solutions we refer the reader to Caffarelli and Salsa [8] . 

Local existence of classical solutions for the classical Stefan problem was shown by Meirmanov 
(see [3S] and references therein) and Hanzawa '33]. In the first approach, the author regularizes 
the problem by adding artificial viscosity to (ll.lb|) and fixes the moving domain by switching to 
so-called von Mises variables. The obtained solutions however, lose derivatives with respect to the 
assumed regularity on the initial data. Similarly, in [33] the author uses Nash-Moser iteration to 
obtain a local-in-time solution, however again with a significant derivative loss with respect to the 
initial data. A local existence result for the one-phase n-dimensional Stefan problem is proved in [33] , 
where the required regularity class for the temperature function is Wp'^ withp>n-|-2. For the two- 
phase Stefan problem a local existence result is presented in |37| in the framework of maximal 
regularity, where the corresponding functional spaces of Sobolev-type require p>n + 3, where n is 
the dimension of the ambient space. 

In a related work, local existence for the two-dimensional two-phase Muskat problem (with both 
different viscosities and different densities) was proved by Cordoba, Cordoba & Gancedo [15] 
and in three dimensions in |16| . Their methods rely on a formulation of the problem as a non-local 
evolution problem for the moving boundary only, whereby the corresponding Taylor sign condition 
is imposed. In a subsequent work [14) . various global existence results were established. An overview 
can be found in [9]. 

As to the Stefan problem with surface tension (also known as the Stefan problem with Gibbs- 
Thomson correction), global weak existence theory (without uniqueness) is developed in [HiSB] 40 . 
In ^27] the authors consider the Stefan problem with small surface tension i.e. ct^ 1 whereby p.lcp 
is replaced by v — an. Local existence of classical solutions is studied in [38]. In [22] the authors 
prove a local existence and uniqueness result for classical solutions under a smallness assumption on 
the initial datum close to flat hypersurfaces . Global existence close to flat hyper-surfaces is proved 
in |32j and close to stationary spheres for the two-phase problem in ,3lj. 

With the Gibbs-Thomson correction, problem (11.11) can account for phenomena such as the phase 
nucleation, undercooling (superheating) and it is also used in modeling crystal growth [37]. It is a 
small-scale model as opposed to the macro-scale classical Stefan problem. In this sense, there is a 
fundamental importance in rigorously understanding the link between the two models. As explained 
in |47j , |46j , one can associate a free energy to the Stefan problem with surface tension defined by 

F,{p,f) = [ pdx + a\r\, 

where p, F are time-independent. Then in the the sense of F-convergence of De Giorgi ^21j . the 
free energy Fo-(p,F) converges to the free energy for the classical Stefan problem, see [47]. This is, 
however, a completely time-independent consideration and does not address the vanishing cr-limit 
of time-dependent solutions to the full non- linear problem (jl.ip . In the context of the water wave 
problem, the vanishing surface tension limit in two and three dimensions has been studied in [2] [3]. 

In the case of the Hcle-Shaw problem (with no surface tension), due to its ill-posed nature, adding 
surface tension can in fact have a singular effect on the behavior of the solution. In particular, even 
in the case of very small surface tension, there is evidence that the vanishing surface tension limit 
does not converge to the solutions of the problem without surface tension [3T1|331I351 130] . Thus, the 
description of surface tension as a stabilizing effect for the system, must be subjected to a rigorous 
verification depending on the problem at hand. In particular, if we think of adding surface tension 
as a singular perturbation of the zero-surface tension problem, the question whether we can establish 
the vanishing surface tension limit becomes a-priori unclear. 
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Turning our attention to the Stefan problem, we can observe that there are two parallel de- 
velopments in the weak existence theory briefly mentioned above. The first one applies to the 
classical Stefan problem and it is motivated by the validity of maximum principle: suitable notions 
of weak/ viscosity solutions have been developed [H [51 |M1 H]. The second development refers 
to the problem with surface tension and state of the art are the existence results in BV-type space 
relying on the gradient flow structure of the problem. From the point of view of the vanishing surface 
tension, it is natural to ask whether the two concepts are compatible in any rigorous mathematical 
manner. The answer is inconclusive for we lack compactness: while the control on the solutions 
from [35l [1] is strong enough to pass to some limit as cr— )-0, it is too weak to guarantee a sharp 
interface in the limit (in other words, to preclude the formation of so-called "mushy" regions [47]). 

We develop a new energy method for the Stefan problem with and without surface tension and 
prove the vanishing surface tension limit. The well-posedness is established in L^-based energy 
spaces determined naturally by the problem. Such a framework is motivated by the analysis of the 
free-surface incompressible Euler equations of Coutand and Shkoller (TB], [TO] . 

Precise statements of our results are given in Theorems ll.ll and ll.2l The estimates that we use are 
non-linear in nature and they fundamentally exploit the intricate energy structure of the problem. In 
particular, no derivative loss occurs with respect to the regularity of the initial data. This framework 
is particularly convenient, as it allows us to rigorously establish the vanishing surface tension limit 
locally-in-time, as formulated in Theorem 1 1.2 1 In this way we link two fundamental models of phase 
transitions that are valid on different spatial scales, thus answering the open question explained 
above. In a forthcoming work, we shall extend our results to the two-phase case providing the 
analogs of Theorems 11.11 and 11.21 

1.6. Methodology and outline of the paper. There are two main ingredients in our approach 
to the Stefan problem. We first introduce the so-called Arbitrary Lagrange-Eulerian variables (from 
now on referred to as ALE-variables) simultaneously imposing a suitable gauge condition on the 
change of coordinates - map 5*. Upon rewriting the problem in the form (jl.9l) . the second main 
input is the extraction of a high-order energy contribution from the moving boundary, a term of the 
general form: 



for k some positive integer. If we were to stick to the original (Eulerian) framework (|1.1[) . the energy 
dissipation law would read 



While formally the same as the basic energy structure for the heat equation, the moving boundary 
contribution exists, but it is only implicitly contained in it. A very important consequence of 
the ALE-formulation of the problem is to make this contribution clear and transparent, in the 
form explicitly given by (jl.2ip . Observe that this energy structure incorporates the Taylor sign 
condition (I1.12p in a very natural way, as it becomes a factor whose sign is crucial to the coercivity 
properties of the quadratic form displayed in (|1.2ip . A further subtlety consists in the discovery 
of another coercive energy term which is defined on the whole domain D, (phase), displayed in the 
fourth line of (|2.33l) . It contains terms of general form 



for a,& as in (j2.33l) . They are intrinsically linked to the problem and carry an information about 
the regularity of the divergence of v] taking a — and b=l the first term above becomes the norm 
of the ALE-divergence of v, as it is easily seen from (|1.9ap : 

The gauge condition (|1.3[) . which relies on a result from [20], allows us to get optimal regularity of 
the temperature function q. This is possible as we have exact (and in general optimal!) control of 
the ff^-Sobolev norms of 4* on the interior, in terms of ff*~'^-^-norm of 4* on the boundary F. This 



(1.21) 





\\d''d^q + d''d^'i'-v\\la.^2 and \\d''d^tqt + d''d^-i'fv\ 
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allows US to prove that the energy S defined in (|2.26p is in fact bounded by the coercive quadratic 
form (the "natural energy" ) F (|2.33p dictated by the Stefan problem. 

Condition is the exact equivalent of the Taylor sign condition, necessary for well-posedness 

of free-surface incompressible Euler equations without surface tension [T5] or the water wave prob- 
lem I48| . If the initial temperature Qq is nonnegative, it is implied by the Hopf 's lemma, at least over 
a short period of time. In a short time regime, we prove a uniform lower bound on A (cf. (I1.13p ). 
thus enabling us to close the estimates. 

In many free-boundary problems, constructing the solution is in general a challenging problem 
despite the (possible) availability of good a-priori estimates. Our main technical idea to make the 
construction as straightforward as possible, is to regularize the problem via horizontal convolution 
by layers as introduced in [TB] in the study of well-posedness of the incompressible Euler equation 
on a moving domain. In addition to that, we also regularize the Stefan condition p = on T{t) by 
modifying it into a Robin-type condition. It is carefully designed to fully honor the energy structure 
of the problem allowing us to prove the basic estimate (Proposition 12 . 5p : 



on its interval of existence {S^ is defined by (|2.26p ). where P is some polynomial. Such a polynomial 
inequality, through a continuity argument leads to uniform-in-K time of existence [0,T] and the 
bound 



Passing to the limit as k— >0, we recover the solution of the Stefan problem (jl.9|) . Our regularization 
is intrinsic to the problem and it does not rely on formulating a sequence of iterated linear problems. 

The second part of this work focuses on the problem of the vanishing surface tension limit. Once 
the well-posedness framework of Theorem ll.ll is set-up, the idea is rather straightforward. Namely, at 
the level of energy, the presence of surface tension simply augments the high-order energy functional 
by a cr-dependent contribution coming from the boundary F, so to obtain (|1.19p . The goal is to prove 
a uniform-in-<T upper bound on on a u-independent time interval [0,T]. This is made possible 
by one fundamental property of S"^ : it distinguishes between two boundary energy contributions of 
the general form 



for suitable a,b€No. Since the error terms are at least of cubic order, we can afford to estimate 
all lower order terms in terms of the a -independent energy term, while the two terms with highest 
number of derivatives get bounded via the a-dependent energy contribution. With uniform estimates 
in hand, we can pass to the limit as ct— ^ 0. 

The plan of the paper is as follows. In Section [2. II we introduce the K-regularized problem and the 
associated high-order energy £k- We then state the energy identities fLemmas 12.21 and 12. 3|) . prove 
that £k is controlled by the "natural energy" (Proposition 12. 4p and finally prove Lemmas 12.21 
and l2.3l In Section [2^ we provide the energy estimates for the error terms. Passage to the vanishing 
surface tension limit is explained in Section |3l In Section |4] we explain how to extend our results to 
the three dimensional setting. 



2.1. A sequence of approximate k problems. We regularize the problem by using the horizontal 
convolution by layers, introduced in 18 in the study of well-posedness of the incompressible Euler 
equation on a moving domain. 

Definition 2 (Horizontal convolution by layers). Let be a C°°(K)-bump function supported in a 
ball of radius k defined through: Pk{x) := j:p{^), where 



£.{t)<C£M + C{t + Vt)P{£^) 



\^^2d''dth\l and a\d'^+'d\h\l 



2. The classical Stefan problem 



(2.22) 
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and constant is such that ^^p{x')dx' = 1. For any given f;:i7— >R we define the horizontal convo- 
lution by layers of g via 

A^g{x^,x'^):= / g{x^ ,x'^)p^{x^ - x')dx' . 



To formulate the regularized problem, we introduce the following quantities: 

'^K(t,x') = {x' ,hi^{t,x')) with hi^{t,x'):=Kf^ki^h{t,x), x' 

and we define ^ on 51 as an extension of its boundary value on F analogously to (|1.3p . In particular, 
■.= deiD'^ ^^\n^(t)\/\n\. Furthermore, 

In analogy to (|1.5p . we introduce 

(2.23) ^2.^ ('9feK -1) . .^.^j-i^ 

Since depends only on t, the Piola identity holds: (51^), fe = for any i = 1,2. The K-regularization 
of the problem (|1.9p takes the form: 



(2.24a) qt-^i{^U,k),j = -v-^u in 

(2.24b) v' + %'iq,k^O in Q; 

(2.24c) q=-K^v-'^al + K'^(3it,x') on F; 

(2.24d) ^'fn« = -wn^ on F; 

(2.24e) t;-A^ = on 917top; 

(2.24f) g(0,•)=g^ *(0,-)=Id, 

where f3{t,x') is defined as 

(2.25) ^(i,^'):^^^a^'=(„.«„2)|^^^^ 

fc=0 

Note that /3{t,x') is a cubic polynomial in t with x'-dependent coefficients. In equation ()2.24dp . 
stands for the outer unit normal with respect to the regularized surface F^, i.e. in the coordinate 
representation 

jdh^-l) 

Tl — 

Note that the corresponding unit tangent to F^ is given via 



Initial values q'^ e C°°(rj) in ^(HM are chosen to converge to go in the energy norm, as k— >0. 

Remark 7 (The regularization). Note that in addition to the horizontal convolution by layers, we 
also altered the Stefan boundary condition q = into its regularized version as displayed in \2.24c\j . 
Two observations are in order. First, this Robin-type regularization is carefully chosen to be com- 
patible with the non-linear structure of the problem. It fully honors the energy structure and the 
smoothing effect is seen in the new energy terms that correspond to n-dependent Sobolev norms of 
ht Ii2.26\) . Second, let us assume for the moment that f3{t,x') = Q. Using I2.24b\l , we would be able 
to rewrite I2.24c^ as 

2 T IK 2|2 I 2 kaI k 2 2 T IK 2|2 , /^f 4\ 

q = n JkI a.l g,2+K q,!^^,- a, = K J„| a,| 9,2 + 0(k ). 

Recalling that ~ 1 for short times, to the leading order in k the above expression shows that q 
is negative on F since (7,2 <0. Therefore, the regularization \2.24c^ with /3 = is consistent with 
the Hopf's lemma which tells us that the negativity of q inside the domain Q implies (7,2 <0 on the 
boundary. 
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Remark 8. We introduce the regularization \2.24c^ to circumvent a technical difficulty of closing 
the energy estimates at the level of highest-in-time differentiated problem. The problem arises from 
the commutation of the horizontal convolution operator in terms of the schematic form: 



j^KK-^tf-^tuT dx', 



where T is a lower order term. Of course, when performing a-priori estimates, such an issue does 
not arise. 

Remark 9 (Compatibility conditions). In order for the solution {q,h) to {2.24-^ to be continuously 
differentiate with respect to t at time t = Q, we must ensure that suitable compatibility conditions on 
the initial temperature q hold. Due to the introduction of k? j3{t,x') term in \2.24r^ , the compatibility 
conditions for the k -regularized problem turn out to be exactly the same as il.l5]) and il.l6]) (holding 
for q'^ instead ofq). 

Theorem 2.1 (Smooth solutions to the K-problem). For k>0 fixed, there is a time depending 
on K, such that there exists a unique solution q = q{K) to \2.24^ on the time interval [OjTk]- Because 
of the compatibility conditions hi. 15]) - {1.11^ , the solution satisfies 

3 

a=0 

Proof. Thanks to (|2.24b ). the boundary condition (j2.24b ) can be written as 

-K^Vq-A'af + q^K^P. 

As such, (|2.24p is a uniformly parabolic equation with C°°-coefficients, with Neumann boundary 
conditions on dVL^op and Robin boundary conditions on F. Standard parabolic theory then provides 
existence and uniqueness of solutions with the stated regularity. □ 

Remark 10. A priori, the time of existence Ti^^O as k— >0. By obtaining n-independent bounds 
on solutions to \2.24^ , we will prove that, in fact, the time of existence is independent of k and given 
byT>0. 

2.2. The higher-order energy function compatible with the k^O asymptotics. The asymp- 
totically consistent higher-order energy function associated to our sequence of regularized K-problems 
is given by 

£,^£^{q,h):= \\d'^d^v\\%^.+ J2 II^^^MliriS 

a+2b<i a-(-2fc<3 

a+2fc<4 a+26<3 
2 12 1 



(2.26) 



h=0 6=0 6=0 6=0 

Remark 11. Thanks to Theorem \2.1[ the function E^, is continuous in time. 

By choosing T^ possibly smaller we assume that for certain (5 > 
(2.27) X{t) > 6 and \dh\l^ < 1/2. 

Since the Taylor sign condition (|1.13l) holds initially, first assumption is guaranteed by the Hopf's 
lemma and the second one follows since dh — initially. 

2.3. The energy identities. In this section we collect the high-order energy identities in two 
lemmas stated below. We use the notation T for those error terms which in an easy straightforward 
way are seen to satisfy the energy bound of the form: 

\T\<P(£.); 

this bound will then always follow from the standard L°° — L'^ — L'^ type estimates. 
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Lemma 2.2. Assume that {q,h) is a solution to the Stefan problem U.9\) (note that a^O). Then 
the following identities hold: 



(2.28) S"'^ latJY ^atj^ 7r 



-1,54, |2 



= / 7^l+ / 7^2+r ; 

(2.29) 



2 



7^3+ / 7^4+r ; 



(2.30) i=i"'^^ zmJy ^atj^^ Jr 
'"7^5 + J Ue+T, 

where TZi, i~l,...6 are error terms given below respectively by {2.4^^ , {2.4^3^ , \2.44^ , \2.45^ , \2.46^ 
and P7^ . 

Lemma 2.3. The following identities hold: 

(2.31) ~i^"''^Jn Jr Jn ^ Jr 
= I Si+ I S2+T; 



EJI / / (-9,2)|M.vI/,.a^P+ / (^Bqu+d^,-vr + ^^ f J-'\dhu\' 

(2.32) i=i_ "'r Jn ^ «rJr 



= / 53+ / S,+T, 
Jn Jr 

where Si, i = l,...4 are error terms given below respectively by \2.52]) . \2.5S]) . \2.54% h2.55\] 

We postpone the proof of Lemmas 12.21 and 12.31 to Section 12. 5[ since we first want to address the 
important question of coercivity of the energy (|2.26p . 

2.4. Coercivity of the energy . By summing the left-hand sides of the identities (|2.28l) - (|2.32p 
from Lemmas 12.21 and 12. 3[ the natural coercive quadratic form that arises as the "energy" takes 
the form 



(2.33) 



o+2b<4 o+2b<3 

W E \J~.'''d'^d%\l.,.^ + ^ E \J'.'~d^dlh,\l^,. 

a+2b<i a+2h<3 
a+2b<i: a+2fc<3 



2 

a+2fc<4; a+2fc<3; 



The mathematical reason for imposing the Taylor sign condition (I1.13P now becomes apparent. In 
order for the second line in the definition of (|2.33p above to make sense we clearly must have 

A(i):= inf (-g,2)(t,x',0)>0, 

x' 
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as it was demanded in (|1.13p . In order to perform the estimates in the next section, it is crucial to 
show that the energy is bounded by P(J-k) for some polynomial P. Such a cocrcivity bound is 
precisely the content of the next proposition. 

Proposition 2.4. The following inequality holds: 

Proof. Due to the Taylor sign assumption (I2.27p . note that the energy contributions from the 
boundary F in the expressions for (|2.26p and for (|2.33p . will satisfy the following bound: 

2 1 
b=0 b=0 

< E \V^2d''d\K^h\l^^.+ \V^2d''dlKht\liLi- 

a+2h<4; a+26<3 

To finish the proof of the proposition we must show that the space-time derivatives of q which appear 
in the energy expression (|2.26l) are also bounded by P(J^k)- Thus it remains to prove the following 
claim: 

Claim. The following estimates hold: 
(a) 

a+2b<i 

(b) 

E W'^'^diqUH^'^n^'^)- 

a+26<4 

Proof of Part (a): It follows from the elliptic theory that \\q\\2l^P{^K) since /S.^^iq — qt + v-w and 
\\(lt\\L=-L^<P{^>i) and \\v\\L^Ll\\w\\LfLl<\\v\\L^Ll\ht\L^Hl<P{^K)- Differentiating (|2.24ap with 
respect to (j = l,2), we obtain /S.^,q,j = [AY^ A^),jq,rnn+{qt + v ■w),j . Furthermore, since q,j = 
^ Tj-v we have 

h^jtW < \\'i',jt\\LrL^\\v\\L^Ll + \\'^,j\\LrL^\\vt\\L^Ll 
< \ht\L^H^\\v\\L^L2+\\D'i'\\L^H^.5\\vt\\L^Ll 

Note that 

<||w||l~hi||w||l~l~ + ||w||l~l2||w,j||l~l~ 

<PiTM\LrH^Jht\L^Hl+J'K\ht\LrHl<PiJ'.)- 

It is easy to see that |j (51™51"),j q,mn li ^ P{J'k)- Thus, by the elliptic theory again, we conclude 

|k||ic.^3<P(J-K). 

Differentiating (|2.24ap with respect to t, we obtain A\s,qt^—{'A^'A"),tq,mn+vtt + {v-'^)t- Again, 
using ||'ytt|||oo^2 ^J'k, the previous estimates and the bound |/it|ioo^2 ^P{J^k), elliptic regularity im- 
plies \\qt\\L^H{<P{J'.)- Furthermore |ktt|iio. ^2 < + '^t • i-lli^.^^ + ll'V;* • < P(J-,). The 
last equality follows from the third line on the definition ()2.26p of J-^ and a simple bound on the 
L^L^-norm of v, which follows from Sobolev embedding. Finally, choose any j,k£ {1,2}. Applying 
d^jd^k to (|2.24ap . we arrive at the elliptic equation 

By the estimates already derived before and the definition (|2.26p , the right-hand side is bounded by 
P(J^k) in L^i^-norm. Thus, by elliptic regularity, we finally conclude ||(?||looj:^4 < P(J'k), concluding 
the proof of part (a) of the lemma. 

Proof of Part (h): The easiest case is 6 = 2. For j' = l,2, we have 

q,jtt = (*K J • w)tt = *K,itt • w -I- 2*K jt • -I- J • wtt . 
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From the above we easily infer 

f\\Dqu\W,< /VtdlLllHlL + ll*«,.tllLll«J§+ll*«,.llLll««llo<^(-^«)- 

Jo Jo 

Since HfZitHo 5 ^ Ikttllo + ll^9ttll-o 5j '^^ conclude the claim when 6 = 2. The rest of the proof relies 
on the estimates established in part (a) of the lemma and the elliptic estimates just like in the proof 
of part (a). □ 



2.5. Proof of Lemmas 12.21 and 12.31 Proof of Lemma \2.S\ Applying the tangential differential 
operator d'^ to the equation (|2.24bp . multiplying it by d'^v^ and integrating over fi, we obtain 

3 

(2.34) {d\^ + d^%1q,k + '^':d\u, d\^) ^.=Y.ci (5'51,^a4-'<z,fc, d\^) , 
where c; — (*). Recalling (|1.8p . we write 

(2.35) B^'A^ = ~^d^^l,,%l + {<9^?lf }, 

where {9'*, 51*^} stands for the lower order commutator defined in (|1.8|) . With this identity, we obtain 

(2.36) =-/ qM'^l~d''^l^r~d'v^N^+ j^'Atd^^l%%kd'v],+T 

Jr ' Jn 

where we used (51|),s = and the identity — ~'A'^q_k to write the last line more concisely. Fur- 
thermore, integrating by parts with respect to x'' 



(2.37) (5lfa\fc,aS^)^,= / A'^B^q.kd^'^ / A'^d'^qd^'-N''- / afa^qaVfe. 

Jn Jr Jn 

Note that the the boundary contribution coming from the fixed boundary d^ltop vanishes due to the 
boundary condition p.24ep . which further reduces to v'^ = on diltop- Summing p.36p and p.37p . 
we obtain 

{d'%^q,k + A^d\k,d\')^,,^^,=- f q,k'Atd'K^^,d\'N^+ I %^,d\d\'-N^ 

(2.38) I- ' > Jr Jr 

Jn 

The first three terms on the right-hand side of (|2.38p will be the source of positive definite quadratic 
contributions to the energy. To extract the quadratic coercive contribution from the first integral 
on the right-hand side of ()2.38p . we simplify it to 

-J q,k'A^d^^l'A';d''v'N' = J {-q,2)Afd^^lAld\'~ J g,i5l2a%;5liaV 

= l^{-q.2 )dH. ■ ^IB^v ■ %l - J^qaBHl^lB^v ■ 

We rewrite the expression B'^v-'A^ and thereby use the boundary condition (|2.24dp : 
B^vAl =B^w%l+B\v~w)-%l 

4 

^B^wAl+d^{{v--w)-'Al)-J2aiB^~'(.v-w)-B'Al 



=0 '=1 



1=1 
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Due to the above identity and recalling ^'^ = A^ A^^*, we obtain 
(2.39) 

4 

>/ r «/ r ^ -j^ >/ r 

The first term on the right-hand side of (|2.39p is rewritten in the following way 
(2.40) 

5^A,* 



= ^dt I {~q,2WA^^>-^l\^^^ I {-q,2tWA,^-'Aiy- / {~q,2)d^A,^-%id^A,^-'Ai 



2|2 



A,9% 



A„ [{-q,2 ) {BHt ■ ~ {~q,2 ) (A,a%t -51^)512 



2^* / (-9>2)|a^A,vl/.5l2|V / A,a% 



+r. 



Ak [(-q,2 ) • ^1)^1] - {^q,2 ) {AjHt ■ ^1)^1 

The second term on the right-hand side of (|2.38p turns into 

%'IB\B\'-N'' = - [ BSB^v-%l^K^ [ B^ht{B\v%l)-Y^aiB'vB^-''Al)+K^ [ p{t,x')B^v%l 

Jr Jr -^r 

= K^ f J-iia^/itp-K^Va, f B^htB'vB^-'^l + K^ f d^p{t,x')d\-%l, 



1=0 



where we used the boundary condition p.24cp in the second equality above (recall ti • 'b^ = w • 'h^ = 
ht). As to the third term on the right-hand side of (|2.38p . note that 

4 4 

'A'^B^v\k = B\%'^v\k ) -J2ciB'%'^B*-'v\k^-B\qt + v 'hi^) -J2r4B'%'^B*~'v\k , 



1=1 



1=1 



where 5lfw*,fe = — div,p^t) = — (qt-f w- V;) by the parabolic equation (|2.24ap . Thus 

af9V,fe(a4g+a%,-«) 

4 

■ I B\qt + -^^fv){d\+d''-^^-v)+y^ci [ d'%'^B''-'v\t{d\+B*^^-v] 
Jn i^-^ Jn 

-\dtl {d\+d^^,,-vy+ f {J2did^-'^.fB'v-B^^.-vt){d\+BH.- 



(2.41) 



1=1 



1=1 J^^ 



Combining ([OS)) . ([O^ . (pillll) and (piiTj) wc obtain the identity with the error terms TZi 

and 7?.2 given by: 

(2.42) 

3 4 4 

n,:=Y,cid'^^B^-'q,kB^v'~{Y,ciB'^'lB^-'v\k+Y.diB^-'^v-B'v-B^^,-vt){B\+BH^ 



1=1 



1=1 



1=1 



7^2 := -A,a% [a, [{~q,2 ) • 51^)51^] - (-g,2 ) (^'A^vft • a^)^^ 
(2.43) (-9,2 )9%K • 5l2a4-'(t; - u;) • ^'^^ -f g,i B*^';,%IB\ ■ %l + ^aza^/z^a't^a^-'a^ 



-LWl3B*v%l. 



i=a 
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Applying the tangential differential operator d'^dt to the equation ()2.24bp . multiplying it by d'^dtv^ 
and integrating over fl, we obtain in a completely analogous fashion identity (12.29^ claimed in 
Lemma 12.21 with error terms TZ^ and TZ^ given by: 

fO AA\ l<m+n<2 l<m+n<2 

a<m+n<2 



(2.45) 



A, [{~q,2 ) {B'dt^t ■ ^1)^1] - (-9,2 ) (d^dtA.^t ■^D'Al 



+ E / {-q,2)d'd't^^-'Aid'~'dr'{v~w)-di'Ai 
i+i'>i •'^ 

0<l+l'<3 



Finally, applying du to the equation (I2.24b[) . multiplying it by duv^ and integrating over ft, the last 
identity (|2.30l) of Lemma 12.21 follows with error terms 72-5 and TZq given by 

(2.46) 7^5 := - (5lf v\k +2^-,t v\kt +2vt ■%t+Vtf'^- ^.tt ■ «t) {qtt + *Ktt • v); 

(2.47) 

7^6 := -A^^u [a„ [(-g,2 • ^D'Al] - (-q,2 ) {dttA^^t ■ 51^)512" 

2 „ 1 

+E«'/ i-q,2)d^'i'.-'AlBt'iv-w)-dl'Al + qaduK^ldttv-%l + K'Y.^uhtdi:vdt'^^ 

1=1 ''^ l'=0 

□ 

Proof of Lemma \2.3l 

Applying the tangential operator B^dt to the equation (I2.24bp . multiplying by B^v^ and integrating 
over we obtain 

{B''dtv\B^v') + {B^dt{nU,k),B\') ^,^0, 

implying 
(2.48) 
1 



-dt / \B^v'\^ + {B''dt^U,k+'A'lB^3tq,k:B' 



^ ci^,jrk{B'd^%'lB'dfq,k.B'v^^ 



1 + 1 = 3, k + k = l 
0<i + fc<4 



RecaUing (jl.Sp . we write 
(2.49) 



Using this decomposition we have 



(2.50) {d'dt'AU,k+^id''dtq,k,d\')^, = - [ %';8'\u,-,%U,kd'v'+ [ ^'ld'dtq,,d'v' +T, 

where the commutator term has been absorbed in the error T. Integrating by parts with respect to 
s and k in the first two integrals on the right-hand side above respectively, we obtain analogously 
to the proof of Lemma 12.21 



(2.51) 



XB''%u,:^U.kB\'+ I %1d'dtq,kB^v' 



+ [ {B^qt + B^%u-vf+ [ {y"diB^-^wB^v + y2eiB^^tv^~^v\,){B^qt + B^'%u-v)+T. 



16 



MAHIR HADZIC AND STEVE SHKOLLER 



Note further that the first term on the right-hand side above can be, similarly to (|2.40p . further 
written as 

J r r 

The second term on the right-hand side of (|2.51D reads, using the boundary condition (|2.24cp 

2 

Jr Jr Jr 

2 

Jr , r. Jr Jr 



Ydi 



1=0 



where the error term T is used to absorb the lower order terms containing the time derivative of J^. 
We also used the regularized boundary condition (|2.24cp in the first equality above. Combining (I2.48P 
- (|2.51l) and the last identity we obtain the identity (j2.3ip with error terms Si and ^2 given by 
(2.52) 

3 3 

Si := E Ci^kJ,k9'd'^A'ld'dfq,k -{^diB^-'w ■ O'v + Y^eid'Afd^-'v^s ) {B^qt+B'^u ■ v) ; 

!+r=3,fc+fc=i 1=1 1=1 

0<! + fc<4 

(2.53) 

52 := -Aj^w ■ [A^[{-q,2){B^w ■ "AD^I] - {-q,2){B^A^w ■ ^1)^1] 

3 2 

+ ^ Q {-q,2 )d^'hj ■ Al {B^-^ iv-w)-B^Al) + q,i AlB^VB^v -Al + K^Y, ciB^BthtB^v ■ B^-^Al 

l=l_ _ 1=0 

-K^d^dtiSd^vAl. 

Applying the tangential operator Bdf to the equation (j2.24bl) . multiplying by BOtv"^ and integrating 
over n we obtain the identity ()2.32p in an analogous way, with error terms ^3 and ^4 given by 



(2.54) 



53 := {v,, -B^^uAt - {BBlA'^}q,k)Bvl + CrnnB''d'tA'lB'-^Bf-'q,k Bvl, 

l<m+n<2 

+ ( J2 d^nB''B^^AlB^"'Bl~\\s - {^tt ■ Bv + B^t ■ vt + ^tBvt )) {Bqtt + B^^u ■ v) ; 



l<m+n<2 



54 := (-9,2 )Bwt ■ Al [{B{v -w)-Al) + (wt -Vt)-BAl + {w~v)- BAI^] 
(2.55) -BBtA^w ■ Al [A^{{-q,2 )AlBdtw ■ Aj) - (-9,2 )B8tA^wAl 

+q,i AlBBt'hv'^BBtv ■ Al + K^Bdtht{B8t {v-Al)~ Bdtv -Al)- K^BBu^Btv ■ A^ 



□ 



2.6. Estimates. The following proposition states the desired energy bound for the classical Stefan 
problem (cr = 0), that will subsequently allow to obtain a uniform-in-K existence interval for our 
family of regularized problems (|2.24p . 

Proposition 2.5. There exists a constant C independent of k and a universal polynomial P such 
that for any td [OjT"] the following bound holds 

£.{t)<C£M + C{t + ^t)P{£,). 

Proof. The proof of the proposition proceeds by systematically estimating error terms in the energy 
identities from Section \2l3[ 
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Estimates for JoJ^Ui. Term: 'ZLiIo In^i^'^td'^''Q,kd'^v' ■ If; = l,wehave 



Jfl 



<|| / ddt{'A'l)\\H^.4d\k\\L^^Ll\\d^V^\\^Ll 
JQ 



For / = 2,3 we have 
rt 



Jn 



< \\d'%'l\\L^H'>.4d^''q,k\\H".^ / l|9V|U. 

Jo 



Term: J2t=i^i /o/o^'^i^^^ ^v^ ,k \d'^q + d'^'if ^ ■ v) . Only cases l — l and I — A deserve special atten- 
tion, while the cases I = 2 and I = 3 are estimates in the usual way. li 1 = 1 we have 



Jo Jn 



+\m'\\Lr 



< 



< / \\d%%.5\\\d\^-'v\kh4d\ + d^'i>.-v\\o 
Jo 

iap-v,fc|U2(||a\llo.5 + ll(a%K-«)llo.5) 

\\q\U4\\dW\0.5 + -Hlo.s) + II fddt{%'l)\ 



IkllL 



Jo 



<^AM\^HlMllHr^+^A>\\LlHl\\D^^ J h\\4.5 + CVt\M\L2Ht' J \\q\\l, 

+cVt\\DHA\L^Ht4D^'^^\\LrH2^-4q\\L-,Hr^<VtP{£.). 

As for the case ? = 4, we estimate as follows: 

t 



Jn 



< 



Jn 



'A'l\\L4\d\°-'[^\k{d\ + dH^-v 



L2 



< J ||i?'*K||2.5|k\fc||wo.r,,oo||a*q + 9>«-i;||o.5<||i?'*K||L-H2.5||q||i=.H3.5 J (klU.S + ||i?'*«||2.5) 

Term: 'Zt=idilQ + IH = 4, we estimate 

^vd^v{d\ + d^^^-v) < f\MU\d\\\o\\dS + d^-^^-v\\o 

Jo 

< Vt\\%,\\LrHyo\\d\ + d^^,-vU^Li\\dM\L-,Li<VtP{£.); 



Jn 



and analogously for / = 3 



/o Jn 

For 1 = 1,2, we have 



d^-d\{d''q+dH,-v) <Vt\M\LrHi-4dS+d^^^-vh^Li\\dMLUi^^tP(^-)- 



Jn 



d^\u-d'^~^v{d^q + d'^'f^-v) < I ||9%||o||a'^"'w||oo||a'^9 + <9'^^'K-w||o 

Jo 

< Vt||'V;|U^~ff2||a4g + a%„-t;|Uj.i2||g||i.^4.5<VtP(£:,). 
The last term on the right-hand side of p.42p is estimated as follows: 

d^^,-vt{d\ + d^-^^-v) <Vi||a"*«||L~i2||94q + a%«-«|Uc.i.||«t||^.^i.5<VtP(f«). 
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Estimates for /q/p7?-2- To estimate the first term on the right-hand side of (|2.43p . we need an 
additional lemma to deal with the commutator term: 

Lemma 2.6 ([19^, Lemma 5.1). For F£ W^'°°(r) and G,9Gei^(r), there is a generic constant C 
independent of k such that 

\A.{FdG) - fAjG\ < C\F\wi.^ir) \G\o, 
where W^'°°{T) denotes the Sobolev space of functions h£L^{T) with weak derivative dh£L°°{T). 



For any i,jg{l,2} set i^ = q,251i G = d^^l and apply Lemma to conclude 



i2|2 



Jo ^ Jo 

< sup \q,2^l:^l\l [ \w\l<Pi£,), 
0<s<t Jo 

where we estimate |u'|^2^3 using the equation (j2.24dl) : 

(2.56) \ht\%H.< f\^'i^Jl + \^h,\Hv■^l))\'<Pi£^) f\\q\\l,<P{E.). 

Jo Jo 

Note that we bounded \v\^ by relating it to its norm over f2 via the trace estimate: 

(2.57) lHi?«3<ii«iii.H3.5<m) rikiii5<m). 



wi.°°(r)l^l3 



Thus 



A^dH 



JT 

< 



A4{~q,2){dHf%l)nl]-{~q,2){d^A^^f^l)^l 



2\ 1/2 



Jr 



Finally, we treat the last term on the right-hand side of (|2.43l) . For 1 < / < 2, we have 



Jt 



< 



- g,2 \l^hi \d'^l\LrHi \d^A, ■ 51^kj^L. VI(|t.|i.ff3 + \w\lihi) < VtP{£.), 



where estimates p.56p and (|2.57l) were used in the last inequality. If 1 — 3, we apply a similar 
estimate, bounding the term — d'^'A, in L^-norm and d'^~\v — w) — d{v — w) via L°° norm and 
Sobolev embedding leading to adesired bound: 



r / {'q,2)d^^.-'Ald{v-w)-d^%l <VtPi£.) 
Jo Jt 



Case Z = 4 is the trickiest error term as four derivatives hit 51^ , thus creating a term that at highest 
order contains five derivatives of '5, being above the energy space regularity allowed for '5. The 
following lemma proves that this term is indeed of lower order, relying on its specific structure and 
integration by parts: 

Lemma 2.7. The following identity holds: 



(2.58) 



- f i-q,2)^d^^,-Ald^^,-d%l+ [ (-g,2)5%«-5l2£;, 
Jt \dT4 Jt 



where E is the error term specified in 112. 59\) below. 

Remark 12. Before we proceed to the proof note that the the function ^! enters the right-hand side 
of the above identity at most with 4 derivatives. This fact enables us to perform the energy estimates. 
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Proof. Note that 

U - W = (w - UJ) • Hk^^k + (w - w) ■ TkTk = (w - • TkTk 

where we used the boundary condition (j2.24dp . We may thus write 

3 

1=1 

where we first used the identity t; • = and in the last hne we used the product rule expansion of 



the the identity = 5 (r^ • with e; the expansion coefficients. Since 



, we have 



1=1 



where we simply used the product rule to expand d {-^-^) and the orthogonality of 9^'^ and 51, in 
the last line. Combining the previous two identities, we may write 

where the error term E is given by 

3 3 

(2.59) E = wT,Y.^id'+'^J^-\\d^,r')-^l + Y,eiWTj'Tj^''^l 



1=1 



1=1 



We thus obtain 



Note that first integral on the right-hand side has a symmetry allowing us to extract a full tangential 
derivative of the level of highest order terms: 



2jr 
1 



\dr.\ 



where we used integration by parts in the second equation. Finally, summing the previous two 
identities we arrive at p.58p . □ 

By standard L°° — LF' — type estimates and Lemma 12.71 we finally arrive at 

(2.60) /* / {-q,-2)d^^.-'Al{v~w)-d''%l <VtP{£.). 

Jo Jr 

Before we estimate the third term on the right-hand side of (12.431) . we first rewrite: 

3 3 



1=0 



1=0 



As a consequence, we have 



(2.61) 







3 , rt 

1=0 



Jr 



10 Jr 
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The first term on the right-hand side above is easily bounded as follows: 



The second term on the right-hand side of (|2.6ip is a sum, and the hardest summand to bound 
is created when l — O. In this case, roughly speaking we bound |9^51^|o by k~^|Ak5'|4 trading one 
tangential derivative on Q^A^A^-D^' for a bound on A^D^' in , at the expense of a factor of . 
Using this observation we obtain 

Jo Jr 

The next-to-last term on the right-hand side of (|2.43l) is again a sum and the hardest term to estimate 
is created again when 1 = 0. We use the same idea as in the previous estimate to obtain 



d^htvd'"At 



Jr 



^\K,^LrLi<ViP{E^). 



Note that we exploited the presence of the K-dependent energy term in our energy f^, using the 
bound K\d'^ht\ 1^-2 <\f£^. In analogous manner, we conclude 



where we note that the commutator term, i.e. the first term on the right-hand side of p.47p deserves 
special attention. Due to the absence of spatial derivatives acting on dtt'^t in 

-A^^u [Ak [(-<Z,2 ) i^ttt ■ Al)%l] - {-q,2 ) {duK^t ■ Al)Al 

we can not apply the commutator bound from Lemma 12.61 in the form stated. Here we crucially 
exploit the ^-dependent term in the energy ■ Note that 
ft 



(2.62) 



A.*, 



A, [(-g,2 ) {^ttt ■ Al) Al\ - i~q,2 ) {duA.^t ■ Al)Al 

<VtP{E^), 



where we gain one power of k in the second line above from the commutator estimate and then 
absorb it into the energy contribution k\'^ ttt\ ■ The last term on the right-hand side of (|2.43p 
contains the /3-contribution from the regularized Dirichlet condition (|2.24cl) . It is easily estimated 
using the Cauchy-Schwarz inequality by a term of the form CtmQ + CtK^\d'^ht\^j^2]^2 which in turn is 
smaller than a constant multiple of tmQ + tEK. Here mo is a constant, which depends only on the 
initial data. 

2.6.1. Estimates for J^Si and /p52. In the first term on the right-hand side of (I2.52p . the harder 
terms to estimate correspond to the cases {l,k) = (2,1) and (Z,fc) = (2,1). If {l,k) = (2,1), then 



dA'^d^dtq,kd'v 



Jn 



< 



\d{dA'^d^v')ddtq,k\ 



Jn 



< 



\\d^A^\U\d''v'\\o\\ddtq,k\\o+ I \\ddtA'^\\oo\\d''v'\\o\\ddtq,k\\o 



<Vi\\d'A4LT^Hl4dML?^Ll\\qt\\L^,H^2 + \\ J dAl{s)d4Hr.4d^v^\\LiLi\\qt\\^m^ 

<Vi\\D'nLrHi4dMLrL4'it\\LiH^^+t\\Dw\\L2^^^^^^ 

<{t + Vt)P{E,). 



Assume now (/,fc) = (2,l), then 

■ d''dtA'ldq,kd\' 



Jn 



<Vt\\q\\LTHl\\dMLrL4M\LlHl<VtP{E,). 
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The second error term is rather straightforward: for any 1 = 2,3 



< 

If / = !, then 



VtML^Hr^0v\\L^Li\\d\t + d''w-v\\L2Ll<VtPi£^). 

<\\DvU^H^^-4wh^H^Jd%+d'w-v\\^2i^.<ViP{£.). 



[ f d'^wdv{d^qt+d' 
Jo Jn 

Similar analysis yields: 



3 , „t 



1=1 



As for the error term ()2.53p . we start by applying Lemma 12.61 to deal with the commutator term. 
For any z,je{l,2} set i^ = q,2!^^5lf , G = d'^w and apply Lemma [Z!6l to obtain 



(2.63) 



f A,[{-q,2){d^w%l)%l]-{-q,2){d^A,w%l)%l] '< f\q,2^l:%l 
Jo ^ ^ Jo 



2 . fvl2|2 



<t sup \q,2^l:^l\l sup \w\l<tP{£,), 

0<s<t 0<s<t 



where, in order to bound |w|L~ij2, we use the equation (|2.24dp analogously to the bound p.56p . 
Upon using Cauchy-Schwarz we get 

f f Aj^w-[A4i-q,2){d^w-%l)%l]-i-q,2)id^A,w-'Al)%l] <VtP{£,). 
Jo Jr 



As for the second term on r.h.s of (|2.53p we obtain 



{~q,2W-\v-w)-d''Al)d''%o-^l < / \q,2\oo{{\v\2 + \w\2W'Al\^)\d^AMo 

Jr 

<\\q\\Lrm--{\\v\\LrHr^ + \w\2)m-id\3 f \d^A^w\o,<ViP{£). 

Jo 

where, the term |w|^oo^2 is bounded by P{£k) for the same reason as in ()2.63p . The last term on 
the right-hand side of (|2.53p is a sum, and the hardest term to bound is created when l = Q. We 
must integrate by parts with respect to the time variable, to obtain 
(2.64) 

t 



K 



10 Jr 
Now observe that 



/ B^htv-B^'^l < ^mo + K\^''ht\LrLl\v\L^L^^^\ I dtd'K\LTLl 



Jo Jr 



^3fyl2| 



< ^^VtK\B%\L2L2<ViP{£.), 

where mg depends only on the initial conditions. As for the remaining three terms on the right-hand 
side of (|2.64l) . they are straightforward to bound using the standard energy estimates. We arrive at 

k"^ f f B^dthtV-B^'Al <K^mo + Vtil + K)Pi£^). 
Jo Jr 

In analogous manner we conclude 



[ Ssdx + [ Sidx' <{t + Vi)P{£.) 
Jn Jr 
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2.7. Proof of Theorem 11.11 The polynomial bound from Proposition 12.51 replaces the typically 
used Gronwall inequality. Using our K-independent estimates together with a continuation argument 
which is detailed in Section 9 of jlT) . there exists a k- independent time T such that 

£n{t)<C£M<C£[0) + l 

for K small enough. Since £ <£k. (recall the definitions (|1.10p of £ and (|2.26p of £k), we obtain the 
uniform bound 

£{q'',h^)<C£{Q) + l, 

where (g'^j/i'*)^ is a family of solutions to the K-regularized problem (|2.24p . < k < 1. Thus, passing 
to the weak limit as k— >0 we obtain a solution on the time interval [0,r] which belongs to the 
space S{T) p.l4p . Since S imbeds compactly into C^'^ the solution is also classical. Note that the 
assumptions (|2.27p remain valid (on a possibly smaller) time interval [0,T], as both |9/i|l^lc^ and 
5 are easily controlled by the energy £. 

Uniqueness. We only present a brief sketch of the uniqueness argument. A simple application 
of the energy method also implies uniqueness of the solution. Assume that {q,h) also solves (|2.24p 
with the corresponding ,v^w. Then the pair (r,p) (q — q,h — h) satisfies the following system of 
equations: 



(2.65a) rt-Al{A'lr,k),j^{A^-^q,)q-iv-i,)-w + v{w-w) in 

(2.65b) {v-vY + A'lr,k+q,kiAf-A'>)^0 in Q; 

(2.65c) r = on T; 

(2.65d) Pt = -r,2 on T; 

(2.65e) dnV — O on dfltop- 



Furthermore, initially (r(0,a;),p(0,a;')) = (0,0). Applying d to the identity (I2.65bp . multiplying by 
{d{v — v)y and integrating over f2, we derive the first identity in analogy to the proof of Lemma l2.2l 
Similarly, applying dt to (|2.65bp , multiplying by (v — vY and integrating, we obtain the second energy 
identity. The natural quadratic form that emerges is equivalent to 

E:^\\d{v-d)\\l,^,+\\v^d\\l^^, + \\r\\l^„,+\\n\\l,^ 

Furthermore, we have an a-priori control of the high-order derivatives of the two solutions, i.e. for 
some M > 0: £{q,h)+£{q,h) < M . From here, we can easily prove the polynomial bound 

E{t)<tP{E{t)), 

which in particular, uses the fact that the initial values for p and r are 0. We infer that E = Q and 
hence the uniqueness follows. □ 

3. The vanishing surface tension limit 

The problem of local-in-time existence for the Stefan problem with surface tension becomes con- 
siderably simpler in comparison to the classical Stefan problem. One important feature of our 
method is that it extends directly to the case with the positive surface tension. Even though there is 
a variety of ways to regularize the problem in the presence of surface tension, for the sake of brevity, 
we shall formulate the K-regularized problem exactly as (|2.24p apart for the condition (|2.24cP , which 
is replaced by 

(3.66) q=-aK,^^^-n\-Al+p on F, 

with (3 defined by (|2.25p . The problem (|2.24l) with p.66p in lieu of (|2.24cp is the K-regularization of 
the Stefan problem with surface tension. 

With cr > 0, we can prove the following energy identities in the same way as Lemmas 12.21 and [231 
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Lemma 3.1. The following identity holds: 

j^^iq,^)it)^ f /{7^l+7^3+7^5+5l+53}+ / / {7^2+7^4+7^6 +52+54} 



Jn Jo Jr 

Jo Jr 



where 

a+2b<A a+2b<3 

andTZi, i~l,...,6, Si, i = l,...,4 are given just like in Lemmas \2.^ and \2.3[ Furthermore, 
(3.67) 

n^^:={-aA^t^),,d^^^-%ld^v-nl + a{~d''A,h-d{-d*A^ht\d^^\-^^ 

+aE«/^-'(^-.).a-'5i^a^A4|M) 

(3.68) +||5^9*A«/^pa,(|a*.rV-i) + aa3aA.«;.51^[aft(|^)-|^] 

+a E au,d'd,A.{^)B'd^;{w-v)-d'~^dr'^l; 
i+i'>i 

(3.69) 

K ( - ,1 ^« • ^Idttv ■ %l + ^\dA^htt\\m.\-''J-^)t - adA^huA,htud{m.\-^) 
(3.70) 

3 — 

S^:=aJ2'^id'-'i^-v)-d''Ald'dtdA^{^)+a E B'^A^w ■^ld^+'d'lAM''''dl-'' {\d^ ,\-') 

1 = 1 a + b<i 

a<3, 6<1 

S^,:=a E a,,,d'-^dr'{w-v)-d^d^%lddudA4^) 

l + l'<2 I '^1 

(3.71) '.''^i 

+CT E &u'aftAK«;-5i^a'+i9,'A./i9i-'a2-' (15*^1-1). 

l + l'<3 
1<1J'<2 

3.1. Estimates. In the following proposition we prove the basic energy estimate in complete analogy 
to Proposition 12.51 In particular, it yields the well-posedness statement under the Taylor sign 
assumption (|1.13p as required in Theorem 11.21 Most importantly, it gives us bounds uniform in a 
that finally allow passage to the limit cr— > 0, necessary to finish the proof of Theorem 1 1.2 1 

Proposition 3.2. There exists a constant C independent of a and a universal polynomial P such 
that for any i G [0,T'^] the following bound holds 

£:{t)<c£:{o)+c{t+vt)p{£:). 

In particular, there exists a time T>0 independent of cr, a constant C* > and the solution {q'^ ,^!'') 
to the Stefan problem with surface tension defined on [0,T] satisfying the bound 

£:'"(g'",^''")<C*, 0<a<l. 
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Proof. In comparison to the estimates for the classical Stefan problem carried over in Section [ 
the only new error terms to estimate are the terms TZ2, Ti-I, TZg, iSf , ^4 given in the statement of 
Lemma 13.11 

Estimating 7?.2 • We start by bounding the first term on the right-hand side of p.67p . 



Jr 



< 



Pi\V^dAM\V^^>.'f\M 



<P{\V^dA^h\L^H^)\^A,'f\L^H^oVt\v\Lm-^<ViP{S^). 

The second and the third term on the right-hand side of (|3.67l) are estimated analogously and rely 
on the standard L°° — L^ — L^ estimates. As for the fourth term on the right-hand side of p.67p . 
note that due to ([2?23| 



Jr 



4 d^Md (^)-^)) 



< 



d^A^h B'^A^h . 



\,<VtP{£^^), 



where the last estimate follows in the standard way: terms with less derivatives are bounded in 
L°°-norm and then further via Sobolev imbedding. In the last term on the right-hand side of p.67p . 
the hardest case to deal with is 1 — 4. Note that 

4 



Jr 



{v~w)-d^nlA4d^A^h\d^^\-^) 

a f [ {v-w)-d^'Aiy^ai,d''~''A^hd''{\d^^\-^)^:I + II. 
Jo Jr 



The more challenging term to estimate is term /. Observe that 



\I\ = a 
+a 

= a 
1 

= r 
1 

= —a 
2 

where 



\v — w 



Jr 



)-d^%lA4d''A^h\d^, 



1-3^ 



<cr 



j-\v-w)-d\dh^-\)d^h^\d^^\-^) 



Jr 



Jr 



J-\v-w)-d\dh^~\)[A4W'A^h\d^^\-^)-W'h^\d^^\~^)) 
J-^~d{\W'h^f){v^-w^)\d^^\-^ 



J~\v-w)-d\dK-\){A4d^AJi\d^^\-^)-d^h^\d^^\-^)) 
commutator term 



commutator term 



Jr 



commutator term := J-i(i;- w) • (A^(96A^/i|9*«|"3) _ a6;^^|^^^|-3)^_ 

In the above we used the parametric representation of 5* in terms of h and integrated by parts. The 
last inequality is rather standard and follows by estimating d{{y^ —w^)\d'^ ^\~^^ in L°° norm and 
further via Sobolev inequality, whereas we also know cr|9^/iK|L°°L2 '^^■'^ ■ To bound the commutator 
term appearing above, we apply the commutator estimate of Lemma l2.6l Term // is easily estimated 
via the standard energy L'^ — l? — l? bound as well as the remaining cases Z = 1,2,3 in estimating 
the fourth term on the right-hand side of p.67p : first integrate by parts so to remove one 9-derivative 
from term and then apply the standard energy estimates. 

Estimating Tl.% . The estimates are completely analogous to the ones for ■ 

Estimating TZq . The first term on the right-hand side of p.69p is estimated analogously to the first 
term on the right-hand side of (|3.67p . Note that 

^ f f \dA,hu\^i\d^,\-^J-^)t <m.\-^J-')t\ooCT f\dA^htt\l<t\dKUdKt\^£-^<tP{£-J, 
^ Jo Jr Jo 

where we use Sobolev inequality and the definition of £^ to infer 

i3/i«ti'oo<ia/iKd?<^(-9,2)i52/i.tp<f:. 
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and similarly 

(3.72) \dK\lo< I {-q,2WK\-'<r^. 



r 

Space-time integrals of the third and fourth term on the right-hand side of (13.69^ are bounded in 
the usual way by P{£^). To bound the last term on the right-hand side of p.69p we distinguish the 
cases I — and ^ = 1. If Z = 1, by Leibniz rule expand 

For the first two terms above integrate by parts to move one d derivative away from d'^A^hu and 
B^Ai^ht- Then use the standard — L'^ — type estimates as well as the bound l^wt 1^=0^2 < 
lkt||L5"ff2.5 to get the desired estimate. For the third term on right-hand side above we have 

°<\d'A,h\^vi\M\d^.r')u\L^^Limv-w)\L^L^^^^^ 

Estimating S2 and SI- The estimates are straightforward and follow the same principle: terms 
with least amount of derivatives are bounded via Sobolev embedding by the cr-independent energy 

Summing up the above estimates we prove the first inequality in the proposition. The existence 
of follows by passing to the limit k— >-0 in the same way as in Subsection 12. 71 The existence 

of (T-independent time T follows from the standard continuity argument and the fact that constant 
C in the proved inequality in the proposition is cr-independent. Since f^(0)— >-£(0) as cr— >-0, the 
last statement of the proposition follows suit. □ 

Proof of Theorem 11.21 Assume that \\{q'' ^'^'^) — [q^ ^'^'^)\\(ji.2 does not converge to as cr->0. 
Then there exists an e>0 and a subsequence (cr„)„gN7 Cn^O as rt— >oo, such that 

(3.73) ||((?'^",vI/-")-(gO,*0)||ci,2>e VneN. 

Since ^'(g'^",^'"'") < C, there exists a subsequence of (g"^" ,^''^'*)n (without loss of generality indexed 
again by {<Jn)) and {q,^) (zS such that 

(g'"",^''"")^(g,*), weakly in 5. 

Note that the injection operator / : 5— > C^'^ is compact. Hence (g'^" ,^''^") — > [q,^] in C^'^ and since 
cr„— >-0 as rt— T'oo, {q,^) solves the classical Stefan problem with initial conditions (g(0),5'(0)) = 
(gojid). From the uniqueness statement of Theorem 11.11 we conclude that (q,^) — {q''^ Thus 
(g'^",^'^-)_^(gO^^O) (ji,2 contradicting ([STTg]) . □ 



4. The three-dimensional case 

In this section, we briefly sketch how to adapt the analysis of the previous sections to prove 
theorems analogous to Theorem 11.11 and 1 1 . 21 in the three dimensional setting. We assume now that 
fl(t) cM^ is an evolving phase in R"^ and the initial and reference domain are given by 

n = T'^x]o,i[, 

where stands for the 2-torus. Initially at i = the bottom boundary 

T = T^x{x^ = 0} 

is the initial reference boundary of the moving domain T{t), whereas the top boundary dCltop — 
X {x^ = 1} is fixed and the temperature p satisfies the homogeneous Neumann boundary condition 
on dfl just like in (|1.2p . We parametrize boundary as a graph over T with the height function h{t,x'), 
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where x' :— {x^,x'^). Via the same ALE change of variables as in (|1.6p . we obtain the fixed boundary 
problem given by (jl.9p . The associated energy is given by: 




|a|+2fc<5 |Q|+2h<4 

(4.74) +1 E \V^2V"d1h\l^^.+ Y IV^V"9f/it|i.^. 

a|+26<4 

q|+26<4; 

In the above definition, a = {ai,a2) is a multi-index of order |a| = ai -|-a2, whereby ai,a2 are non- 
negative integers. Symbol V refers to differentiation in tangential directions, i.e. ^"■=d'^ld"i. The 
three-dimensional Taylor sign condition for a function q reads: 

(4.75) A3^(t):= inf (-g,3)(i,a;',0) >0. 
The following theorem holds: 

Theorem 4.1. Let the initial conditions {qo,hQ) be such that £^^{qt),ho) <cxi and let qo satisfy 
the Taylor sign condition \4-75^ . Then the three-dimensional one-phase classical Stefan problem is 
locally-in-time well-posed, i.e. there is aT>0 such that there exists a unique solution {q,h) with the 
initial data {qo,ho) on the time interval [0,T]. In addition it satisfies the bound: 

£:3^(g,/i)<2f3^(go,/^o). 

Furthermore, let (gQ,^Q)(j>o be a given family of well-prepared initial conditions in the sense of Defi,- 
nition\^ Assume that it satisfies the Taylor sign condition \4.75^ and the corresponding compatibility 
conditions. By (q"^ ,h")cr>o we denote the associated family of solutions to the problem lll.9\) . There 
exists a a-independent time T>0 and a constant C depending only on (qo^ho) such that 

£^^'''{q'',h''){T)<C a>0. 

for all a>0. As a consequence, sequence {q" ,h'^) converges to the unique solution {q,h) of the 
classical Stefan problem \1.9\) with a = in C^'^ norm. 

Remark 13. Note that the definition of £^^ contains time derivatives. Thus, to make sense out 
of the assumption £^^(qQ,hQ) <oo, we express the time derivatives dtqo and dth^ in terms of the 
spatial derivatives as explained in Remark\^ 
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